A field-oriented chain of dipolar particles in elongational flow 
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We study the behavior of an isolated field-oriented chain of dipolar particles in elongational fluid 
flow. Our main goal is to emphasize the effect of dipolar interactions on the chain's contribution to 
the pressure tensor and to the viscosities of a dilute suspension of these linear aggregates. In our 
model, despite the overall rigid appearance of the chain at rest, the constituent beads may move 
slightly relative to one another, conferring a certain degree of flexibility to the chain. This flexibility 
is quantified in terms of a dimensionless parameter, A -1 , comparing thermal and dipolar energies. 
We perform an expansion in A" 1 , and obtain the first correction to the rigid chain contribution to 
the Kramers' pressure tensor for different flow geometries. The interplay of the elongational flow 
field and the field-induced chain orientation gives rise to a rich variety of scenarios. We compute 
the elongational, shear, and rotational viscosities in some representative situations. 
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I. INTRODUCTION 

Colloidal suspensions play an important role in many 
natural phenomena as well as in various industrial pro- 
cesses. The stability of these suspensions against aggre- 
gation of their constituent particles is of essential impor- 
tance for their behavioiu. As long as interparticle in- 
teractions can be neglected, one does not observe aggre- 
gation and the characterization of the system reduces to 
understanding how the physical properties of the fluid, — 
e.g. its viscosity — are modified due to the presence of the 
particlesEni On the contrary, the phenomenology of the 
suspension may change dramatically if interactions be- 
tween particles become relevant. 

In the late 1930's, WinslowO first noted the curious be- 
havior of a suspension of dielectric particles in oil when 
subjected to an electric field. He reported the forma- 
tion of linear chains of particles aligned with the electric 
field, and how the effective viscosity of the suspension 
could change by orders of magnitude by simply modify- 
ing the applied field. An analogous field-induced behav- 
ior is exhibited by magnetorheological fluids, e.g., by a 
suspension of magnetizable superparamagnetic particles 
in a nonmagnetic fluidlll, or by a suspension of aqnmagne- 
tizable spheres in a ferrofluid (magnetic holes pEJ. Since 
then, several experimental and theoretical efforts, as well 
as computer simulations, have been devoted to study the 
different aspects of the complex-behavior of colloidal sus- 
pensions of dipolar particleaBilj. 

The change in the rheological properties of dipolar sus- 
pensions upon the action of an external field, is in part 
due to the aggregation of the colloidal particles, which 
form clusters of macroscopic size. These are usually lin- 
ear chains — quasi rigid rods — oriented along the direc- 



tion of the applied field, although for high enough con- 
centrations, of-dipolar particles more complex structures 
may ariseE3~E3. The overall spatial arrangement of the 
aggregates is very effective in hindering the fluid flow, 
conferring the suspension a solid-like texture. 

Our purpose in this paper is to ascertain the role played 
by strong, but finite, dipolar interactions in the linear 
rheology of a suspension containing either an isolated 
field-induced aggregate or, equivalently, a dilute concen- 
tration of such clusters. After introducing the dimension- 
less parameter, A, comparing dipolar and thermal ener- 
gies, we obtain the first correction in powers of A -1 to 
the rigid rod (A -1 — > 0) contribution to the Kramers' 
pressure tensor for different flow geometries. In particu- 
lar, we analyze the influence that those interactions exert 
on the value of some of the transport coefficients of the 
system. We have organized the paper in the following 
way: Section || contains a description of the system and 
the conditions under study. The end-to-end vector of the 
chain is introduced and its mean value, as well as the 
mean square end-to-end distance, are obtained. In Sec- 
tion [II we analyze the contribution of a chain, under 



the action of an elongational flow, to the pressure tensor 
of the suspension, and to some of the viscosities charac- 
terizing the system for different elongational flow fields 
of interest. The conclusions are summed up in the last 
section. 



II. DESCRIPTION OF THE SYSTEM 

The structure of a colloidal particle aggregate may be 
either fixed or deformable, depending upon the nature 
of the aggregation process and the type of interparticle 
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bonds established. In Ref. [u], the issue of the growth 
of similar electric field-induced chains was experimen- 
tally addressed. There, the authors reported that, in 
the presence of an intense field, the colloidal dipoles or- 
ganize themselves into linear chains. These clusters ex- 
hibit a certain degree of flexibility, limited of course by 
the strength of the dipolar interaction among the beads. 
Within one chain, they observed fluctuations in the sep- 
aration between sphere surfaces (5 — 10% of a diameter), 
as well as in the angle between the relative position vec- 
tor of consecutive spheres and the applied field. Hence, 
we will assume that the beads forming the chain are in 
close contact, but still able to slightly move in relation 
to one another due to thermal fluctuations and to the 
action of an elongational flow field. 

The dimensionless parameters describing the condi- 
tions of the system under study are essentially: A = 
m 2 /(d 3 kB T), comparing dipolar and thermal energies, 
and the Langevin parameter /i = mH/kgT, where m 
and d are the dipolar moment and the diameter of each 
sphere, H is the external (magnetic or electric) field 
strength, kg is the.-JBoltzmann constant, and T the ab- 
solute temperatures. In our analysis A is assumed to 
be large and fx — ► oo, as in the experimental realizations 
in which the irreversible formation of these sort of quasi 
rigid linear aggregates is observed. 

Despite the long-range contributions to the energy of a 
chain containing a fixed number of particles, TV, the am- 
plitude of the vibrations and oscillations of the spheres 
are almost entirely taken into account by just considering 
interactions between nearest-neighborsE3l23. Long range 
effects are subleading corrections to the short range at- 
tractive contributions and can be ignored. As a matter 
of fact, the most important contributions correspond to 
relative distances among the particles of the order of one 
diameter, and relative orientations with an azimuthal an- 
gle, defined with respect to the external field direction, 
close to = 0. 

Thus we will assume that the potential energy can be 
expressed as a sum of nearest neighbors terms 
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is the interaction energy of the dipoles i and j, with 
r ij = r ij^ij the vector giving their relative position. Its 
modulus is the distance between the sphere centers. 
Since the most important contributions come from dis- 
tances ~ d and relative orientations 6ij ~ 0, each 
term in the, sum of Eq. (Jj]) can be approximately written 
as followsEJ 




FIG. 1. Chain of magnetic particles in an elongational flow. 
All the magnetic moments m point towards the direction of 
the external field applied. The vector R is the end-to-end vec- 
tor giving the relative position of the centers of the spheres 
on the extremes. 

In the last equation, we have introduced the relative posi- 
tion vector between two consecutive spheres in the chain 
qi = Ri+i — Ri through its components in spherical co- 
ordinates qi = (qi,9i,(fii). The position vector Ri deter- 
mines the center of the i-th. sphere. Moreover, we have 
performed an expansion around the minimum of the en- 
ergy in the variables 6i and £j, which, in turn, is defined 
by qi = d(l + and such that < & -C 1. 

Multipolc and multibody contributions are not consid- 
ered because they are expected to be very small. Like- 
wise, we disregard London- Van der Waals forces, and any 
possible steric interactions due to the surfactant mono- 
layer with which the colloidal particles may be usually 
coated. 



Mean size of the chain 

To characterize the mean size of a chain, containing 
TV dipoles, under the action of an elongational flow and 
thermal fluctuations, let us consider the end-to-end vec- 
tor R giving the relative position of the centers of the 
spheres located on the tips of the chain (see Fig. [l]) 
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In order to compute the average of this vector, or equiva- 
lently of the vector q^, we need the probability distribu- 
tion ip({Ri}) of the spheres' position. This probability is 
the solution of the diffusion equation describing the coa- 
servation of system points in the configurational spaceEa, 
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where fij^ is the relative translational mobility tensor, 
V ma g is the dipole-dipole magnetic potential, and f3 is 
the elongational rate tensor (J3 = f3 T ) corresponding to 
the stationary homogeneous external flow, v Q = (3 ■ r. 
As an approximation, we assume that the spheres move 
through the solvent without disturbing the velocity field 
{free draining), so that 
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where, for the sake of simplicity, we consider that all the 
spheres in the chain have the same radius a = d/2, and 
I is the unit tensor. 

The next step is to write the diffusion equation in terms 
of the relative position vectors, qi. To this purpose we 
introduce the transformations 



Qi — £ BikRk, 
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B ik = 5. 
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from which one infers the relation d/dRi 
J2k Bki d/dqk', and the matrix 

f 2 if i = j 

A v = £ B ' kB i k = { - 1 if 1 = j ± 1 • 

k I otherwise 



(8) 



Taking into account Eqs. the diffusion equation 

can be alternatively written as 



dip _ k B T 
dt 6nr]a 

d 



d 



£ Q q . "E^fe ( f) n , 



&l> , dV mag /k B T 



dq k 



dq k 



(9) 



The stationary solution of this equation with a homo- 
geneous potential flow field is 
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where the symbol : indicates a double contraction of 
indexes, and Cy is the Kramers matrix defined in 
Rcf. |25|. Up to linear order in the elongational rate (3 — 
linear dynamics — , we expand the exponential factor in 
Eq. ( |l0| ) . The stationary average of q k is then given by 

N-l 



qk) = / ( H d( li) Qk ipst - de z 
J i=i 
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where we have introduced the translational diffusion coef- 
ficient for a single particle D = ksT/ (6irrja). This result 



is obtained after decomposing the integral into different 
parts, and taking into account the following relations^: 
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In the Appendix we indicate in more detail some of the 
calculations involved in the computation of a similar av- 
erage which appears in the following section. 

Therefore, up to first order in A -1 and in the elonga- 
tional flow rate, the average end-to-end vector R defined 
in Eq. ([|) turns out to be 

(R) ~ (N - l)de z + N(N 2 - IJ-^La- 1 ^ . e z - f3 zz e z ), 
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where, once again, use has been made of the relations 
(|l2l). At equilibrium (/3 = 0), this average obviously re- 
duces to (R)eq — (N — l)de z — the length of a rigid chain. 

Following the same procedure we can compute the 
mean square end-to-end distance 
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Up to first order in (3, and in the small parameter A -1 of 
our expansion, we do not find any contribution from the 
elongational rate. Nevertheless, the equilibrium value is 
given by 



{R 2 ) eq ~(N-l) z d z + 
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Thermal fluctuations effectively stretch a quasi rigid 
field-oriented chain — using the root mean square end- 
to-end distance as a measure of its length — by a global 
amount [(JV - l)/(3A)] 1 / 2 d. On the other hand, ther- 
mal effects, would be obviously negligible for a rigid chain 
(A iV)E§. The stretching factor per bond is 



[(R% q - (R)U^ 
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= [3X(N - l)p 1/2 d. 



(16) 



Thus, for large but finite values of XN, i.e. (AA) _1 ^ 0, 
the mean separation between consecutive spheres in- 
creases by a factor proportional to (AN)^ 1 / 2 . This result 
( |l6| ) is in qualitative and quantitative agreement with the 
observations reported in Ref. |l4| for similar electric field- 
induced chains. In these experiments, A ~ (20 — 30) and 
the mean number of particles in a chain at the late stages 
of the aggregation process is typically N ~ (5 — 15) which, 
according to our result (|l6|), yield a mean near- neighbors 
separation of about (3 — 6)% of d. This value is compa- 
rable to their observations, in which the sphere surfaces 
appear to be separated by (5 — 10)% of a sphere diameter 
due to thermal fluctuations. 
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III. CONTRIBUTION OF THE CHAIN TO THE 
PRESSURE TENSOR 

Besides the contribution of the solvent to the total 
pressure tensor of the suspension, there is another con- 
tribution coming from the direct interaction of the par- 
ticles constituting the chains. The latter can be ob- 
tained from the rhcological equation of state proposed 
by Kramers 
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where V is the volume of the system. 

Once we know the stationary solution of the diffusion 
equation (|lo|), the average appearing on the right hand 
side of equation ([l?]) can be expressed as an integral of 
the form: 
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As we are interested in finding the Newtonian viscos- 
ity tensor, we again expand the exponential factor in 
Eq. (|h]) up to first order in (3. To evaluate the remain- 
ing integral, it is convenient to decompose it into different 
parts as indicated in the Appendix; thereby, in the New- 
tonian domain, the average in Eq. (n8[) reads 



9a- 



k B T\l+—Y, dkP ■ e z e z 



dqk 

d 2 

-77^ A 1 C kk f3 • (7 - e z e z ] 
6D 



(19) 



Finally, after using Eqs. (fl^), we obtain the contribution 
of one chain to the pressure tensor (|17|) 
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The term proportional to (XN)^ 1 constitutes the first 
correction to the rigid chain limit (A -1 — » 0). Depending 
upon the structure of the flow rate, this term represents 
either an increase or a decrease in the value of the dif- 
ferent components of the pressure tensor as a function 
of A, i.e. when the chain becomes more flexible. Note 
that, as we are not taking into account hydrodynamic 
interactions between the particles {free draining approx- 
imation), the rigid chain- limiting pressure tensor grows 
simply as the third power of the chain length. Hydro- 
dynamic interactions are responsible for the logarithmic 
term that should also.-appear when dealing with a long 
straight line of spherescJ, but their effects are less impor- 
tant for shorter chains as the ones considered here. 



Moreover, from this expression we can obtain the con- 
tribution of the chain to the viscosity tensor, rj p , which 
we can identify by comparing with the relationship 
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This tensor is thus 
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where S is the fourth-rank tensor Siju = (SikSji + 
diidjk)/^, symmetric in any of these index transforma- 
tions (i «-> j, k <-> I, (ij) «-> (kl)). 

For the conditions under consideration, the chains are 
always oriented in the direction of the field, but smoothly 
vibrate and oscillate around this orientation. In this case, 
the symmetries of the fluid flow essentially determine the 
characteristics of both the pressure and the viscosity ten- 
sors. For instance, as the chains are oriented along the 
z-axis — parallel to the external field — if the flow field rate 
P is diagonal, so will be the pressure tensor. On the other 
hand, for a non-diagonal flow rate, we have, in general, 
both symmetric and antisymmetric contributions. Let us 
illustrate the structure of the pressure tensor by consid- 
ering some representative cases of interest: 

i) Flow through a pore. In this case, if the system has 
rotational symmetry around the z-axis, the velocity field 
is given by 
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where (3 is the elongational rate. The pressure tensor has 
then the following form 
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Under these conditions the chain is subjected to tensile 
and compressive forces giving rise to the so-called elon- 
gational viscosities defined from the differences between 
two of the diagonal components of the pressure tensor, 
i.e. normal pressure differences. As an example, 
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This quantity essentially grows with the third power of 
the rigid chain length (Nd) 3 , and for high but finite val- 
ues of A. When A -1 ^ 0, the chain becomes slightly more 
flexible, its effective length increases giving rise to a small 
increment in the viscosity of the system. 

Alternatively, a fluid flow with rotational symmetry 
around the x (or y) axis leads to the following structure 
of the pressure tensor: 
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Combinations of any two diagonal components of this 
tensor amount to different elongational viscosity coeffi- 
cients. Moreover, in any of the above mentioned geome- 
tries, we have another viscosity coefficient coming from 
the trace of the pressure tensor. 

ii) Planar elongational flow. This flow can be gener- 
ated by four rotating cylinders. If we locate the cylin- 
ders such that the elongational tensor is again a diagonal 
matrix, in view of Eq. (^p|), the pressure tensor is also 
diagonal and, of course, symmetric. For example, if 



= -fix 



= (3z, 
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we find a similar phenomenology as the one described 
above. On the other hand, in a reference frame rotated 
45° with respect to the previous one, the elongational 
tensor is 
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resulting in 
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which has both symmetric and antisymmetric contribu- 
tions. Related to these parts, we find not only elonga- 
tional viscosities but also a shear viscosity, rj, and a rota- 
tional viscosity, r] r . These last quantities can be identified 
from the non-diagonal components of the symmetric and 
antisymmetric parts, TL%^ and Tlxz > respectively. They 
turn out to be 
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Thus, for A -1 — > 0, both shear and rotational viscosi- 
ties behave like those of a dilute suspension of cylinders 
and grow with the third power of the rigid chain length, 
(Nd) 3 . For high but finite values of A, however, the 
shear viscosity increases, whereas the rotational viscos- 
ity decreases. The reason why the latter occurs is that in 
this particular geometry, the elongational flow not only 
stretches the chain but also tends to reorient it. This 
rotation is opposed by the presence of the external field; 
in the case of a rigid chain (A -1 — > 0), the magnetic field 
impedes more effectively the reorientation and, as a con- 
sequence, the rotational viscosity attains its maximum 



value. Nonetheless, the rotational viscosity is smaller for 
a finite value of the parameter A" 1 , or, in other words, 
for a slightly more flexible chain. In Fig. || we plot the 
relative variation of the shear viscosity as a function of 
the number of particles in the chain N, for three given 
values of A -1 = 0.01,0.05,0.1 within the experimental 
range. If we define the shear viscosity in the A -1 — > 
limit as r? , the relative correction (r]—r) )/r) — (3AiV) _1 
grows linearly with A , and becomes less important for 
longer chains. The relative decrease of the rotational vis- 
cosity \r) r — r^\/rj^ shows the same behavior depicted in 

Fig. a 



I 




FIG. 2. Relative correction to the shear viscosity of a chain 
as a function of the number of dipolar particles TV, for fixed 
values of A" 1 = 0.01, 0.05, 0.1. 



IV. CONCLUSIONS 

In summary, we have computed the contribution of a 
linear chain of dipolar particles to the Kramers' pressure 
tensor and the viscosities of a dilute suspension of such 
clusters. The field-induced chain is oriented along a pre- 
ferred direction — parallel to the external field — but, at 
the same time, is exposed to both thermal effects and 
the influence of an external elongational fluid flow. 

As was pointed out in previous experimental studies, 
these linear clusters exhibit a certain degree of flexibility 
despite its rod-like appearance. At equilibrium, thermal 
fluctuations are responsible for average separations be- 
tween the spheres surfaces of the order of (5 — 10)% of a 
diameter, as well as for slight angular deviations, which 
amount to an effectively larger mean square end-to-end 
length of these field-oriented chains. The flexibility of a 
chain is obviously limited by the strength of the dipo- 
lar interactions among its constituent colloidal particles. 
The dimensionless parameter A, comparing dipolar and 
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thermal energies, thus becomes an essential ingredient in 
our analysis. 

The rigid chain-like rheologic behavior, which corre- 
sponds to the limit A -1 — > 0, is modified considerably 
as the chains become more flexible, i.e., for large but fi- 
nite values of A. Firstly, we have characterized the mean 
size of the chain obtaining the average end-to-end vector 
and the mean square end-to-end distance as a function 
of the parameter A -1 . Secondly, we have computed the 
contribution of the chain to the pressure tensor of the sys- 
tem from the rheological equation of state proposed by 
Kramers. The structure of the pressure tensor has been 
explicitly illustrated for different elongational flows. Fur- 
thermore, by comparing our result to the standard linear 
law relating the pressure tensor and the elongational flow 
rate, we have calculated the viscosity tensor of the sys- 
tem. 

Some of the most common viscosity coefficients char- 
acterizing a colloidal suspension of this kind have also 
been obtained. We have discussed the role played by the 
parameter A in the behavior of these quantities for two 
practical situations of interest: a flow through a pore and 
a planar elongational flow. In the first one, the symmetry 
of the contribution of the chains to the pressure tensor of 
the suspension only gives rise to elongational viscosities, 
defined from the differences between any pair of diago- 
nal components of the pressure tensor. As the flow field 



effectively stretches the chain by a global amount propor- 
tional to A -1 , the elongational viscosity also increases by 
a similar factor beyond its rigid limit value. In the sec- 
ond situation described, the chain's contribution to the 
pressure tensor is no longer symmetric; in this case, we 
have computed the so-called shear and rotational viscosi- 
ties. The former increases with A -1 as well, whereas the 
latter decreases, i.e., the reorientation of the chains is 
slightly favored for large but finite values of this param- 
eter, yielding a rotational viscosity smaller than in the 
case of a rigid rod. 

Our results, valid up to the orders specified in the anal- 
ysis, should provide a good description of the most fre- 
quent experimental conditions and, consequently, could 
be directly contrasted with experiments. 
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APPENDIX: 



In this appendix, we indicate in some detail some of the steps involved in the evaluation of Eq. (|18|). First, we split 
up the integral into different parts, a)-e), which are 
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Let us explicitly compute some of these contributions. For instance, 
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The other contributions are obtained in a similar fashion. Adding up the different parts and considering the relations 
( |l2| ) , we can rewrite the average in Eq. ( |l8| ) as 
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Finally, after the appropriate normalization of the probability density, and up to first order in /3, we obtain 



<?fc- 



dVmac, 

dqk 



= k B T< 

d 2 . 
+ 12D' 

x/l- 



d 2 (3 Zi 
2AD 



■N(N 2 - 1) 



d 2 



(3 : (I-e z e z ) 



N 2 



1 



d 2 (i zz 
24D 



6 

d 2 
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N(N ~ 1] ~ 72D X {N ~ 1] 13 : {I ~ 



which can be simplified to the final form (|l9|). 
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